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We consider a restricted three body problem, where two interacted particles are located in two 
dimensional (2D) plane and interact with the third one located in the parallel spatially separated 
■ plane. The system of such type can be formed in the semiconductor coupled quantum wells, where 

' the electrons (or holes) and direct excitons spatially separated in different parallel neighboring 

^ I quantum wells that are sufficiently close to interact and form negative X~ or positive indirect 

Q , trions. It is shown that at large interwell separations, when the interwell separation much greater 

• than the exciton Bohr radius, this problem can be solved analytically using the cluster approach. 

' Analytical results for the energy spectrum and the wave functions of the spatially indirect trion are 

obtained, their dependence on the interwell separations is analyzed and a conditional probability 
04 ' distribution is calculated. The formation of 2D Wigner crystal of trions at the low densities is 

predicted. It is shown that the critical density of the formation of the trion Wigner crystal is 
sufficiently greater than the critical density of the electron Wigner crystal. 
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PACS numbers: 71.35.Pq, 71.35.-y, 73.21.Fg 



INTRODUCTION 



The linear and nonlinear properties of semiconductor heterostructures are often governed by excitons that are 
defined as the bounded state of the electron-hole pair. This is especially important for semiconductors with wide- 
C ' bandgap, where the excitonic binding energy is comparable to room temperature. Positively {X'^) and negatively 
, charged excitons, called trions in semiconductor structures, have been the subject of many experimental and 

(~| theoretical studies in the last years. Positively charges excitons are formed by one electron and two holes, and 
O ', negatively charged excitons are formed by two electrons and one hole. The theoretical proof of the stability of trions 
in bulk semiconductors was presented by Lampert [ll. The experimental observation of trions in semiconductor 



quantum wells (QWs) was achieved in CdTe/CdZnTe [2| and in GaAs/AlGaAs [3 
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1 The many-body formalism including Feynman diagram technique was developed to study the collective properties 
J> ' of three-dimensional (3D) trions in the bulk semiconductors 0. The 3D and X~ confined in a semiconductor 
C , cylindrical quantum dot were investigated using a variational procedure within the effective mass approximation Q . 
• The 3D trions in bulk semiconductors 's'l and two-dimensional trions Q were studied theoretically. According to these 
studies, the 2D trions have binding energies that are larger than the trions in the corresponding bulk materials. The 2D 
excitons and 2D X^ m QWs in the presence of 2D electron gas were observed experimentally by optical measurements 
of an excitonic recombination line on the photoluminescence spectra from the QW [lo| . The experiments devoted 
to the observation of 2D X^ in QWs in magnetic field by studying the magneto-optical absorption spectra were 
performed in Ref. Qj- Besides, the 2D X' in QWs in magnetic field were observed as a results of the analysis of the 
excitonic recombination line on the electroluminescence spectra from the Q W 11 ill . The experiments devoted to the 
observation of magneto-optical spectra of 2D X~ in QWs were performed 0, [l^l- -R-^f. Q the 2D X~ are studied 
in coupled quantum wells, where electrons and holes are spatially separated in different QWs. 
; ^ ' The problem of finding the wave functions and eigenenergies of X~ and X^ formed by three particles is the tree-body 
I problem. There are different approaches were applied to solve three-body problem: the numerical diagonalization 
of Hamiltonian for the eigenvalue problem, the variational approach, the Hartree-Fock method, the approach of 
Faddeev equations, the method of the hyperspherical functions (HF). The three-electron quantum dot in a harmonic 
confinement potential was studied using the HF method [l^ [l3| and Faddeev equations in the configuration space [l5| . 
The hyperspherical functions method was applied to solve the Schrodinger equation for 2D X+ and X~ [l6|. The 
three-body problem for three electrons in a quantum dot in a magnetic field was solved in the framework of the 
hyperspherical functions method [l^ ■ The variational method was applied to solve the eigenvalues and eigenfucntions 
problem for 3D X+ and X~ in cylindrical quantum dot Q. A variational calculation of the ground-state energy of 



2D neutral excitons, X+, and X~ in a single QW was performed in Ref. [18|. A variational calculation of the lower 
singlet and triplet states of positively charged 2D trions, confined to a single quantum well and in the presence of a 
perpendicular magnetic field, was performed |19|]. The ground-state energy of 2D X+ in coupled QWs (CQWs) formed 
by spatially separated two 2D heavy holes in one QW and one 2D electron in neighboring QW was calculated using 
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the variational method within the infinite-hole-mass approximation ^2^. The 2D X~ in the presence of perpendicular 
magnetic field were studied numerically by the diagonalization of the Hamiltonian in Refs. pTI.[2^. The ground state 
energy of the 2D in the presence of 2D electron gas was calculated using the diagonalization of the eigenvalue 
problem in the random phase approximation ^3^. In Ref. '23'] the screening of the Coulomb interaction and the 
Pauli exclusion principle were considered. A first-principle path integral Monte Carlo study of the binding energy 
of excitons, 2D positively and negatively charged trions and biexcitons bound to single-island interface defects in 
2D GaAs/AlGaAs quantum wells was presented in Ref. [23]. The preliminary results for trions in the system of 
semiconductor coupled quantum wells were reported at the 21 European Conference on few-body problem in physics 

. . . 

In all theoretical studies mentioned above the Schrodinger equation for the trion was solved numerically using the 
different methods. In this Paper we propose the type of the trion which wave functions and eigenenergies can be 
obtained analytically for some conditions. We consider two parallel CQWs. In the case 1 there are electrons in one 
QW and the electrons and holes in the other QW. In the case 2 there are holes in one QW and the electrons and 
holes in the other QW. In the case 1 we consider X~ formed by one electron in one QW and an electron and a hole in 
the other QW. In the case 2 we consider X^ formed by one hole in one QW and an electron and a hole in the other 
QW. The eigenfunction and eigenvalue problem for these trions is the restricted three-body problem, since the motion 
of the electron (or hole) is restricted in the plane of one QW, and the motion of the electron-hole pair is restricted 
in the plane of the other QW. We show that for relatively large distance between these CQWs the eigenvalue and 
eigenfucntion problem for this trions can be solved analytically. 

The dilute system of trions at the densities, when the average distance between the trions is much larger than the 
radius of each trion the system of trions can be treated as the dilute system of the fermions with the pair Coulomb 
repulsion, which can form non-electronic Wigner crystal. The 3D electron Wigner crystal was described in Ref. [2^ 
and the 2D electron Wigner crystal was described in Ref. [l^ . The 2D electron Wigner crystal formed by electrons on 
films of helium was studied in Refs. [1^ H^. The spectrum of collective excitations for the 2D electron Wigner crystal 
formed by electrons on a helium film was calculated [1^. The ground state of a 2D electron gas including formation 
of the 2D Wigner crystal at low densities was calculated using Monte-Carlo method "ss]. Wigner crystallization of 
electrons in 2D quantum dots was analyzed using path integral Monte Carlo approach [34] . We predict the formation 
of 2D Wigner crystal of trions at the low densities. In our consideration we neglect the influence of the disorder in 
the quantum wells. We will show that the critical density of the formation of the trion Wigner crystal is sufficiently 
larger than the critical density of the electron Wigner crystal, since the mass of trion is larger than the mass of the 
electron. 

The Paper is organized in the following way. In Sec.|lT]we introduce the model Hamiltonian for X~ and X^ in the 
CQWs reducing the restricted three-body problem to the 2D three-body problem for the exciton and the projection 
of the electron or hole on the plane of the excitonic QW. In Sec. IIIII we review the solution of the Schrodinger for 2D 
exciton in a single QW. In Sec. IIVI we find the wave function of the relative motion of the exciton and electron (hole) 
and obtain the corresponding energy spectrum. The energy spectrum for the trion, it dependence on the interwell 
separation and calculation of the conditional probability distribution based on the trion wave functions are presented 
in Sec. |Vl The formation of the trion Wigner crystal in the dilute trion gas is discussed in Sec. IVIl Finally, our 
conclusions follow in Sec. IVIIl 



II. THE TOTAL HAMILTONIAN AND WAVE FUNCTION FOR AN ISOLATED TRION 



We consider two parallel coupled quantum wells separated by a spatial barrier of thickness D. There are excitons 
in one QW and electrons or holes in the other QW. One method for obtaining a system of excitons is to create 
excitons in one region of the QW using laser pumping. 2D electrons or holes can be obtained in the QW by doping 
semiconductors that tends to create a constant Fermi energy, producing an excess of carriers, either electrons (donors) 
or holes (acceptors) . These excess carriers then give the material either electrons in the conduction bands or holes in 
the valence band. 

Let ri be the coordinates of particles in laboratory frame with effective masses of the electron mg and the hole mh- 
Let's introduce Jacobi coordinates which are an appropriate set for the three-body problem: 



77 = (r2 - ri) 

/ meYi+mhY2\ 
Ye{h) = ■ 1) 
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FIG. 1: The restricted three-body problem in the CQWs: an electron in one QW and an electron-hole pair in the other QW. 

where mg(^-) is the electron (hole) mass, the coordinate vector t] in 2D space describes of the relative motion of the 
electron and hole in the exciton with the reduce mass fj, = rriemh/ {me + Trih) in the excitonic QW, the coordinate 
vector Ve(h) from the center of mass of the exciton to the electron (hole) in the other QW describes the motion of the 
particles with reduce mass m = (me -I- mh)me(i^^ / (rrie + rrih + mg^^-j), R is the vector of the center of mass of trion 
and M = me -\- mh + me{h) is the mass of a negatively or positively charged trion. 

The total Hamiltonian for an exciton in one QW interacted with an electron (hole) in the other QW is given by 



H = -- 



2m 



e{h) 



-^e(h) + Hex + Ve{h)-ex , 



(2) 



where Ve(h)-ex is the potential energy of the interaction between the electron (hole) in one QW and the exciton in an 
excitonic QW, and Hex is the exciton Hamiltonian. 

Since the motion of the exciton is restricted in the excitonic QW and the motion of the electron (hole) is restricted 
in the electronic QW, we replace the coordinate vector of the electron (hole) 'fCe{h) by its projection r on plane of the 

1/2 

excitonic QW using the relation Te^h) = (^^ + -D^) . Thus, we reduced the restricted 3D three-body problem to the 
2D three-body problem on the plane of the excitonic QW as shown in Fig. [1] 

The Hamiltonian of the electron and the hole forming the exciton in the excitonic QW has the following form 



2me 



2mh 



e 



(3) 



e is the electron charge. The problem of the in-plane motion of interacting electron and 
me and m/j, respectively, forming the exciton in the excitonic QW can be reduced to that 

moving in a Coulomb potential and motion of the center of mass 



where and are the Laplacian operators with respect to the components of the vectors ri and Y2 respectively, 
e is the dielectric constant, 
hole with effective masses 

of one particle with reduced mass /i and energy Eex , 
of the exciton. 

Let us assume, that the distance between QWs D is greater than the Bohr radius as of the exciton formed by 
the bounded state of the electron and the hole in the excitonic QW: D ^ as, where as = ti^e/{2mexe^)- Under 
this assumption the potential energy Ve[h)-ex of the pair interactions between the exciton and electron placed in two 
parallel quantum wells which is the sum of the Coulomb electron-electron repulsion and electron-hole attraction for the 
X~ trion or the hole-hole repulsion and electron-hole attraction for the trion can be reduced to the polarization 
energy. Follow Ref. [35| the polarization energy for the interaction of the isolated 2D exciton and an isolated electron 
can be written in an approximation in the form 



ye{h)-ex{r) ^4^^^ - (^2+^2)2 ' W 

where o; = |i ^-2^ and r is the distance between the center of mass of the exciton and the projection of the electron 
(hole) on the excitonic plane as shown in Fig. [1] In other words, r is the projection of the coordinate re(h) defined by 
Eq. ID) onto the excitonic plane. 



4 



Based on Eqs. ([2]) - Q after the separation of the motion of the center of mass of the trion the eigenstates of the 
trion in coordinates ([IJ are described by the fohowing Schrodinger equation 

-J^A, - ^ - |lAr - j-,r^-^) nv,r) = E ^{r,,r) , (5) 



2n ery 2m {r^ + D^) 

where and Ap are the Laplacian operators with respect to the components of the vectors -q, and r, respectively. 
Eq. ([5]) ahows the separation of variables and reduces to the Schrodinger equation that describes the eigenstates of 
a 2D exciton, and the Schrodinger equation that describes the relative motion of the exciton and electron or exciton 
and hole in the 2D excitonic plane. 

Seeking the solution of (O in the form ^'(r], r) —'^ ^xiv)^ e{h){'^) and using the separation of variables we obtain 



2/i e-q 



^v- — ] = -^exiri), (6) 



^2 a \ 

Therefore, the energy spectrum of the trion is 

E^E,, + E,(^n)- (8) 

Let us consider ^ and ([7]) separately. 



III. THE WAVE FUNCTION AND ENERGY SPECTRUM OF THE 2D EXCITON 



The eigenstates of the electron and the hole forming a 2D exciton in the excitonic QW are described by the 
Schrodinger ([6]), that represents a two-body problem which is like the two-dimensional hydrogenic problem. Below 
we follow the approach of Ref. [36j for finding the wavefunctions and eigenenergies of the 2D hydrogen atom in polar 
coordinates for the arbitrary values of two masses. The Schrodinger equation (j6]) that describes the relative 2D motion 
of an electron and hole, in polar coordinates has the form 



2/^ ydrf rj drj 



£7] 



(9) 



After the standard separation of the variables r] and (p we obtain the radial Schrodinger equation and the solution 
for the 2D exciton wave function in terms of associated Laguerre polynomials is presented in (36j . and the exciton 
spectrum is given by 



Er, 



Ep 



lie 



2e2(n- 1/2)2^' 



(10) 



where n = 1,2,3,... are the quantum numbers. We can see that in the ground state at n = 1 the binding energy 



of the exciton is Ei — 2 fie / (e h ), and the characteristic radius of the exciton corresponding to the wave function 
^^exiv) at n = 1, Z = is given by gb — Tt^e/ (2fie'^). Interestingly enough that by replacing (n — 1/2) by n retrieves 
relevant formulas of the 3D case, particularly quantization conditions ([TU]). 



IV. THE WAVE FUNCTION AND ENERGY OF THE RELATIVE MOTION OF THE EXCITON AND 

ELECTRON AND EXCITON AND HOLE 

The eigenstates of the electron (hole) in the field of the exciton are described by the Schrodinger ([7]) , which represents 
a motion of the projection of the electron (hole) on the plane of the excitonic QW. Using the polar coordinates r and 
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(j) and substituting e{h){f,<t') = f ^''^We(ft,)(r, 0) into Eq. ([7]), we obtain the 2D Schrodinger equation of the relative 
motion of an exciton and an electron (hole) 



'2m 



Qj.2 



1 

4^,2 



r 



1 92 

) +^<h)-ex{r) 



Ue(h) (r, 0) = EexUe{h) (r, (j)) 



(11) 



where Ve[h)^ex{''') is given by Eq. Q. The 2D motion of the projection of the electron (hole) in the plane of the 
excitonic QW, can be treated as that of a single particle with reduced mass m and energy -E'e(^), moving in a potential 
{r'^+D'^y^ ■ Assuming r ^ D, the Taylor expansion of Ve[h)-ex{'>') from Eq. ^ up to the linear term results in 



, - a 1 a 2 



(12) 



£)4 ' 2£)6 

Using potential (|12l) in ([TT|) we obtain the Schrodinger's equation that describes the two-dimensional motion of 

a single electron (hole), that corresponds to the 2D harmonic oscillator with the harmonic parabolic potential 7^, 

where 7 = a/D^, when the energy is counted from the energy level -E^^^-i = —a/D'^. 

Using the standard separation of the variables r and (f) in and follow Refs. (37l.[38l| we obtain the radial Schrodinger 
equation and the solution for the eigenfunctions for the projection of the electron (hole) on the plane of the excitonic 

.1/2 



QW can be written in terms of associated Laguerre polynomials. The characteristic radius a = (h/ (2.^7717)) of 



the wave function is directly proportional to _D^/^ (a - 



-1/2 



D3/2y rpj^g corresponding energy spectrum is given by 



E 



N L 



i?e(.)=-|f + (2iV+l + |L|);.(^) 



1/2 



(13) 



where N — min(n, n'), L — n — n' , n, n' — 0, 1,2,3,... are the quantum numbers. 

As it seen from Eq. (fT3|) . the first term is proportional to D~^, and the second term is proportional to D~^. 
Therefore, for the some values of the separation distance D between the QW the energy becomes positive. That means 
that at the interwell separations greater than some critical Dcr, D > Dcr the energy Epf l is positive, Ej^ ^ > 0, 
and at Z? < Dcr, the energy En l is negative, En l < 0. For example, at the lowest quantum state N — L — 
as it follows from Eq. (|13l) the ground state energy for the electron (hole) in the parabolic well of the exciton in the 
excitonic QW is given by 



Eoo — — 



1/2 



Using Eq. (|14|) . and the value for a we obtain the critical interwell separation Dcr- 

Dcr 



32 ) ?iei/2 



(14) 



(15) 



If we consider the trions in GaAs/AlGaAs CQWs, then we have rric — 0.07mo, nih — O.lSwo, e — 13, where mo is the 
free electron mass. Using Eq. (|15l) for the GaAs/AlGaAs CQWs we obtain for X~ the critical interwell separation 
Dcr = 4.4 nm, and for we find Dcr = 5.6 nm. However, let us mention that our approach is valid only for 
D 3> Dcr, because the condition of the Taylor expansion of the electron (hole)-exciton interaction potential energy is 
D > UB- 



V. DISCUSSION OF THE RESULTS FOR THE SINGLE TRION 

In our cluster approach, when the trion is the system of interacting exciton and electron or hole, the wavefunction 
of the trion \l/(77, r) is a product of the wavefunctions of the exciton and relative motion of exciton and electron (hole) 
and, therefore, can be expressed in terms of associated Laguerre polynomials. As it follows from (|8]) the energy 
spectrum of the trion is defined by Eqs. ([TU]) and and is given by 



En 



N L 



21 e^a% 



fj,e 



32 e£|4 2e2(n- l/2)2?i^ 



+ {2N + l+\L\)h 



/21 e'^a% 
yMemD^ 



1/2 



(16) 



The energy spectrum of the trion corresponds to the discrete levels corresponding to the bounded states of the 
electron (or hole) and exciton. Therefore, the trion wavefunction represents the bounded state with the characteristic 



6 



Energy Spectrum of X* 

n.1, N=0, L.O 

n.1, N=0, L=1 

— - n.1, N=t, L.0 
—■- n.1, N=1, L=1 



Energy spectrum of X ' 

n.1, N.0, L.O 

— - n.1, N.0, L.I 
— - n./, N.1, L.0 
n=1, N.1, L.I 



FIG. 2: The energy spectrum for X and X as a function of the interwell separation D for different quantum states. 



radius at large D corresponding to the radius of the trion at = a. Since a is directly proportional to _D^/^, which is 
the highest degree of D compare to the other characteristic radiuses, thus the radius of the trion is proportional to 
1)3/2, rpj^g calculated energy spectrum as a function of the interwell separation D for different quantum numbers for 
and X" are shown in Fig. [5] 



TABLE I: The energy spectrum of X'^ and X trions for interwell separation D — 12 nm for different quantum states. 
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Energy, meV 
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-12.7 


1 





1 


-11.6 


1 





2 


-10.8 


1 





2 


-9.0 


1 





3 


-8.8 


1 





3 


-6.5 


1 





4 


-6.8 


1 





4 


-3.9 


1 





5 


-4.8 


1 





5 


-1.3 


1 





6 


-2.9 










1 





7 


-0.9 










1 


1 





-10.8 


1 


1 





-9.8 


1 


1 


1 


-8.8 


1 


1 


1 


-6.5 


1 


1 


2 


-6.8 


1 


1 


2 


-3.9 


1 


1 


3 


-4.8 


1 


1 


3 


-1.3 


1 


1 


4 


-2.8 










1 


1 


5 


-0.9 










1 


2 





-6.8 


1 


2 





-3.9 


1 


2 


1 


-4.8 


1 


2 


1 


-1.3 


1 


2 


2 


-2.9 










1 


2 


3 


-0.9 
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Eq. (|16l) shows the dependence of the energy spectrum of a trion on the interwell separation D. Typical thickness of 
the barrier between two parallel QWs D varies from one experiment to the other experiment. The thickness of barrier 
separation for production of indirect excitons was 4 nm (39l. |40|. However, for interwell distances used in the drag 
experiments in Ref. j4i| are D = 17.5 nm and D = 22.5 nm, and the interwell separations used in the drag experiment 
in Ref. [42| are 20 nm and 30 nm. Therefore, we can conclude that it is possible to experimentally construct the QWs 
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with thickness of the barrier from 4 nm to 30 nm. In our calculations we vary the interwell separation within these 
limits (from 10 nm to 20 nm). Analysis of the energy spectrum of the trion, for example for D = 12 nm, shows that 
X+ exists when the exciton is in the ground state with n — 1 and the hole is in states: iV = 0, L = 0, 1, . . . , 7 or 
the exciton is in state n — 2 and the hole is in the ground state. X~ exists when the exciton is in the ground state 
and the electron is in states: iV = 0, L = 0, 1, . . . , 5 or the exciton is in the state n = 2 and the hole is in the ground 
state. Let us once again mention that there is the critical interwell separations when for D > Dcr the the ground 
state energy of the trion is negative, and if we consider the trions in GaAs/AlGaAs coupled QWs, then we obtain for 
X~ the critical interwell separation Dcr = 4.4 nm, and for we find Dcr = 5.6 nm. However, our approach is valid 
only for D ^ Dcr, because the condition of the Taylor expansion of the electron (hole)-exciton interaction potential 
energy is D ^ as- 

It is interesting to mention that from Eq. ([16]) follows that, for example, for L = the maximum value of iV = 3 
for distance D = 12 nm and N = 7 for _D = 15 nm when the exciton is in ground state and the energy of the trion 
is still negative. The energies spectrum of and X^ trions for the interwell separation D ^ 12 nm for different 
quantum numbers are presented in the Table I. 

Analysis of Table I shows that the X^ trion does not exist in the states {n, N, L) = 
(1,0,6), (1,0,7), (1,1,4), (1,1,5)), (1,2,2), (1,2,3), (1,3,0), (1,3,1), while the the X+ trion exists in all 
these states. It is clear from Eq. (1161) that the energy spectrum depends on the trion reduce mass m. The reduced 
mass for the X+ trion is 0.89mo and it is greater than one for the X^ trion, which is 0.53mo, because the effective 
hole mass m/j is greater than the effective electron mass rric, and, therefore, the second positive term for the trion 
energy in Eq. (fTB)) is greater for X~ than for X~^ , since the second term is proportional to m~^/^. Therefore, for the 
states listed above, the binding energy for X^ is positive, and the bound state for X~ does not exist, while for the 
same quantum states the binding energy for X^ is negative, and the bound state for X~^ exists. 

The structure of the wave function of the trion which is a product of the wavefunctions of the exciton and relative 
motion of exciton and electron (hole) can be tested by calculation of the conditional probability distribution (CPD) 
which is the pair correlation function P(r,ro) that expresses the probability of finding an electron or hole at the 
position r given when the exciton is located at the position rg. Therefore, the CPD gives us information on when 
exciton located at rg sees an electron (hole) at r and defined as (30l - [3^ 



In Eq. (jl7p 5* is the wave function of the trion, which is a product of the wavefunctions of the exciton and relative 
motion of exciton and electron (hole). We use the accurate analytical wave functions for the trion to compute the 
CPDs. In Fig.|3]the CPDs are displayed for the state n = l,l = 0, N = 2, L = for the different interwell separations 
that vary from D — 10 nm to I? = 18 nm. At the interwell separation D = 10 nm the confinement potential is strong 
and the contribution of the single-particle energies to the total energy of the trion is larger. The hole is confined in a 
rather compact region so that the CPD does not show clear two peaks. On the contrary, for D — lA nm the effect of 
the confinement becomes weaker. The size of the system grows and we see clearly well separated hump structure in 
Fig. [3}; and Fig. [3ji. For the comparison in Fig. |4] the CPDs for X^ trion is presented when the hole is in the excited 
state with iV = 3 for different interwell separations. The CPDs for four specific states i.e., the n = 1, I = 0, N = 1, 
L ^ 0, the n ^ 1, I = 0, N = 1, L ^ 1, the n ^ 1, I = 0, N ^ 1, L ^ 2 and n = 1, I = 0, N = 1, L = 3 when 
the interwell separation D = 14 nm in conjunction with identification of regularities related to the quantum number 
L are shown in Fig. [5j By examining the CPD's associated with selected states one can make a conclusion on the 
dependence of the CPD on the quantum number L. 



It is well known that the dilute system of the electrons can form the 2D Wigner crystal 33] . Based on the existence 
of the trions in the CQWs we can propose that trions form of 2D Wigner crystal at the low densities alike to the 2D 
Wigner crystal formed by the dilute system of the electrons. The dilute system of trions, when the average distance 
between the trions is much larger than the radius of each trion, can be treated as the dilute system of the fermions 
with the pair Coulomb repulsion, that undergo a phase transition and crystallize by forming the Wigner crystal. To 
minimize the potential energy, the trions experiencing Coulomb repulsion form a triangular lattice in 2D system. In 
other words, this is an example when Wigner crystal phase is occurring in non-electronic system at low density. 

The Wigner crystal formed by the composite particles was studied in several recent publications. Ref. [43| was 
devoted to Wigner crytallization of the dipole cold atoms, while Wigner crystallization of dipole indirect excitons in 




(17) 



VI. TRION WIGNER CRYSTAL 
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FIG. 4: The CPDs for trion for the state n = l, l = Q, N = 'i, L = Q and different inrewell separations: (a) D = 12 nm; 
(b) D = f4 nm; (c) D = 16 nm; (d) D = 18 nm. 
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QWs was studied in Refs. [4J-|46|. Wigner crystallization of the dipole cold atoms in magnetic field was considered 



in Ref. [45[- The Hamiltonian of the trions in CQWs will be written alike to the Hamiltonian of composite dipole 
particles in the papers listed above, when the energy of the composite particles is counted relative to the binding 
energy of the composite dipole particle. The total Hamiltonian Htot per unit area of the interacting system of N 
trions per unit area is similar to that for the electrons [33| : 

Htot = K + U . (18) 

In Eq. (fTS]) K is the kinetic energy of the centers mass of the trions given by 

where S is the area of the system, M_ = rrie + rrih + me is the mass of the trion X~ and — me + m^ + m^ is the 
mass of the trion , and U is the potential energy of the trion-trion Coulomb repulsion per unit area: 

l<i<j<N 

In Eq. (PH)) Ri and Rj are coordinates of the centers of the mass of the trions. Let us mention that in the dilute 
system of trions when the average distance between the trions is much larger than the radius of each trion, we can 
neglect the interaction between excitons of different trions and exciton in one trion and electron (or hole) in the 
other neighbored trion, which are proportional 1/r^ (see Ref. [i^l) and l/r'* (see Eq. ([4])), respectively, compare to 
the Coulomb repulsion between the electrons (or holes) in the different trions. In this approximation we replace the 
Coulomb repulsion between the electrons (or holes) in different trions by the Coulomb repulsion between the centers 
of mass of different trions. Therefore, we can treat the the dilute system of the charged spin- 1/2 fermions interacting 
with the pair Coulomb repulsion formed by the trions the same way as the dilute system of electrons. 

According to Ref. [il], the qualitative criterion of the stability of the Wigner crystal is the condition, when the 
potential energy dominates the kinetic energy: (U) > {K), where {U) is the average potential energy per unit area, 
and {K) is the average value of the kinetic energies of the centers of mass of trions per unit area. This allows to 
determine the density at which the trion system becomes a Wigner crystal. Let us mention that the averaging of 
U and K is calculated as the averaging by the many-particle antisymmetric wave function of the centers of mass of 
trions which is similar to the antisymmetric many-electron wave function [ssj . 
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Let's show that at zero temperature trion Wigner crystal exist at higher densities than electron (or hole) Wigner 
crystal. Estimating tj and K for the trion system analogously to the electron system [s^, we conclude that the 2D 
trion Wigner crystal is stable when the dimensionless density parameter rs > 37 at T = 0. The so-called Wigner-Seitz 
radius is defined as = a/ag, where a = (tt/o)"^/^ is the average distance between the centers of mass of trions, p is 
the 2D density of trions, and ao = ?i^e/(M+(_)e^). The last condition enables us to determine the density at which 
the trion gas becomes a Wigner crystal. For the negative trions X~ this condition corresponds to p < 4.14 x 10^^ m^^ 
and for the positive trions X+ p < 6.75 x 10^"^ m~^. For the electron Wigner crystal we substitute ao = ?i^e/(mee^), 
which results in p < 2.41 x 10^^ m~^. For the hole Wigner crystal we substitute oq — Tt^ e / irrihe^) , which results 
in p < 1.11 X lO^'^ m^^. Therefore, the trion Wigner crystal can be formed at the sufficiently higher densities than 
electron (or hole) Wigner crystal, because the mass of trion is greater than the mass of the electron or hole. 

Above we considered the trion gas at zero temperature. At finite temperatures the defects will always destroy the 
long-range translational order in the 2D Wigner crystals. We assume that a 2D Wigner crystal in the ground state 
will have no defects. Also let us mentioned that in practice, it is difficult to experimentally realize a Wigner crystal 
because quantum mechanical fluctuations overpower the Coulomb repulsion and quickly cause disorder. A low dense 
system is needed. 

VII. CONCLUSIONS 

In this Paper we have reduced the three-body restricted problem for a spatially separated exciton and electron 
(hole) located in the CQWs to the 2D three-body problem of the exciton and the projection of the electron (hole) on 
the plane of the exciton QW. In the limit of large spatial separation between the electron (hole) and exciton QWs we 
have obtained the analytical expression for the wave functions of the trions and X+ with the spatially separated 
electron (hole) and exciton. Analytical results for the energy spectrum of the trion formed by the exciton and electron 
(hole) separated in different quantum wells are obtained, its dependence on the interwell separations is analyzed and 
a conditional probability distribution is calculated. The differences in the binding energies for and X^ due to the 
difference between electron and hole masses is analyzed. The 2D Wigner crystallization of the trions in the CQWs 
is predicted. Due to the fact that mass of the trion is greater than the mass of electron, the critical density of the 
formation of the trion Wigner crystal at zero temperature is sufficiently greater than the critical density of the electron 
Wigner crystal. 
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